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Abstract. We apply some recent developments of Baldoni-Beck-Cochet- 
Vergne IBBC V05I on vector partition function, to Kostant's and Steinberg's 
formulae, for classical Lie algebras A r , B r , C r , D r . We therefore get efficient 
Maple programs that compute for these Lie algebras: the multiplicity of a 
weight in an irreducible finite-dimensional representation; the decomposition 
coefficients of the tensor product of two irreducible finite-dimensional repre- 
sentations. These programs can also calculate associated Ehrhart quasipoly- 
nomials. 

Nous appliquons des resultats recents de Baldoni-Beck-Cochet-Vergne B BCV05| 
sur la fonction de partition vectorielle, aux formules de Kostant et de Stein- 
berg, dans le cas des algebres de Lie classiques A r , B r , Cr, D r . Ceci donne 
lieu a des programmes Maple efficaces qui calculent pour ces algebres de 
Lie : la multiplicity d'un poids dans une representation irreductible de di- 
mension finie ; les coefficients de decomposition du produit tensoriel de deux 
representations irreductibles de dimension finie. Ces programmes permettent 
egalement d'evaluer les quasipolynomes d'Ehrhart assocics. 

1. Introduction 

In this note, we are interested in the two following computational problems 
for classical Lie algebras A r , B r , C r , D r : 

• The multiplicity of the weight /z in the representation V(A) of highest 
weight A. 

• Littlewood-Richardson coefficients, that is the multiplicity of the 
representation V{y) in the tensor product of representations of highest 
weights A and /x. 

Softwares HE (from van Leeuwen et al. vL94j) and GAP |GAP| h and Maple 
packages coxeter/weyl (from Stembridgc S95 ), use Freudcnthal's and Klimyk's 
formulae, and work for any semi-simple Lie algebra (not only for classical Lie 
algebras). Unfortunately, these formulae are really sensitive to the size of coeffi- 
cients of weights. Moreover, they do not lead to the computation of associated 
quasipolynomials (A, /x) i— > and (A, /j,, v) i— > c^" . 

Here the approach to these two problems is through vector partition function, 
that is the function computing the number of ways one can decompose a vector as 
a linear combination with nonnegative integral coefficients of a fixed set of vectors. 
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For example the number p{x) of ways of counting x euros with coins, that is 

p{x) = (J{n E Z s + ; x = 200ni + 100n 2 + 50n 3 + 20n 4 + 10n 5 + 5n 6 + 2n 7 + n s }, 

can be seen as the partition of the 1-dimcnsional vector (x) with respects to the set 
{(200), (100), (50), (20), (10), (5), (2), (1)} of 1- dimensional vectors. In the case of 
the decomposition with respects to the set of positive roots of a simple Lie algebra, 
we speak of Kostant partition function. 

Recall that any d-dimensional rational convex polytope can be written as the 
set P($, a) of nonnegative solutions x = (xi) £ R N of an equation Yli=i Xi(f>i = a, 
for a matrix $ with columns <pi E V and a € V (d = N — r). It follows that 
evaluating the vector partition is equivalent to computing the number of integral 
points in a rational convex polytope. 

The vector partition function arises in many areas of mathematics: representa- 
tion theory, flows in networks, magic squares, statistics, crystal bases of quantum 
groups. Its complexity is polynomial in the size of input when the dimension of 
the polytope is fixed, and NP-hard if it can vary |B94l IB971IBP99] , 

There are several approaches to the vector partition problem. For example 
Barvinok's decomposition algorithm B94], recently implemented by the LattE 
team DHTY031 \Q, works for general sets of vectors. Beck-Pixton BP03 also 
created an algorithm dedicated to the vector set arising from the Birkhoff polytope, 
counting the number of semi-magic squares. 

In this note, we use recent results of Baldoni-Bcck-Cochet-Vergne |BBCV05] 
to obtain a fast algorithm for Kostant partition function via inverse Laplace for- 
mula. These results involve DeConcini-Procesi's maximal nested sets (or in short 
MNSs DCP04 ) and iterated residues of rational functions computed by formal 
power series development. 

We combine resulting procedures with Kostant 's and Steinberg's formulae giv- 
ing and c^ v in terms of vector partition function. We then obtain a Maple pro- 
gram computing for classical Lie algebras (A r , B ri C r , D r ), the multiplicity of a 
weight in an irreducible finite-dimensional representation, as well as decomposition 
coefficients of the tensor product of two irreducible finite-dimensional representa- 
tions. To the best of our knowledge, they are also the only ones able to compute 
associated piecewise-defined quasipolynomials (A,//) i— » and (A, v) i— ► c^. 

These programs (available at [C]) are specially designed for large parameters 
of weights. Indeed although only written in Maple they can perform examples with 
weights with 5 digits coordinates, far beyond classical softwares written in C++. 
We also stress that our programs are absolutely clear, easy to use and require no 
installation of exotic package or program. Retro-compatibility has been checked 
downto Maple Vr5. They are fully commented, so that a curious user can figure 
out their internal mechanisms. 

However, certain other softwares and packages are not limited by the rank of 
the algebra like our programs. For example computation of non-trivial examples 
in Lie algebras of rank 10 is possible with the software HE, whereas our programs 
are efficient up to rank 5-7. These facts make our programs complementary to 
traditional softwares. 

Remark that Kostant 's and Steinberg's formulae have already been imple- 
mented once in the case of A r C03 . This previous program relies on results 
of Baldoni-Vergne |BV01| implemented by Baldoni-DeLoera-Vergne |BdLV03] . 
computing Kostant partition function only in the case of A r . Tools were special 
permutations and again iterated residues of rational fraction. 
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A new technique for Littlewood- Richardson coefficients has been recently de- 
signed by DeLoera-McAllister jDM05| . For A r , they wrote an algorithm using 
hive polytopes |KT99| . For B r , C r , D r , they implemented Berenstein-Zelevinsky 
polytopes BZOl]. They can also evaluate stretched Littlewood- Richardson coeffi- 
cients c t £ v tjl . These two methods consist in computing a tensor product coefficient 
as the number of lattice points in just one specific convex rational polytope. How- 
ever our programs based on multidimensional residues are faster, and can reach 
examples not available by their method. 

This paper is organized as follows. Section |21 recalls representation theory 
problems we are interested in and links them with algebraic combinatorics. Sec- 
tion describes more precisely rational convex polytopes and formulae counting 
their integral points. Section 01 introduces maximal nested sets and formulae that 
were used in our programs. Finally in Section[5]we perform tests of our programs. 

2. Representation theory and convex polytopes 

Let us fix the notations once and for all. Let g be a semi-simple Lie algebra 
of rank r. Choose a Cartan subalgebra t of g and denote by L C f* the weight 
lattice. 

Let A + be a positive roots system. The root lattice is defined as Z[A + ]. Let 
C(A + ) be the cone spanned by linear combinations with nonnegative coefficients 
of positive roots. The Weyl group of g for t is denoted by W. 

There exist only four simple Lie algebras A ri B r , C r , D r of rank r, called 
classical Lie algebras of rank r Bou68 , and determined by their positive roots 
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The character of a representation V of g is ch.(V) = X^eL dim(V Al )e M . Recall 
that the irreducible finite-dimensional representation of g of highest weight A is 
denoted by V(X). Hence the weight multiplicity is defined as dim(T/(A) At ) for 
any weight (i such that A — fi is in the root lattice. Multiplicities are called 
Kostka numbers when = ^4,.= s[ r+1 (C). 

On the other hand, multiplicities of representations V(y) in the tensor prod- 
uct V(A) ®V{[i) are called Littlewood- Richards on coefficients (or Clebsch-Gordan 
coefficients). Here v is a dominant weight such that X + fi — v is in the root lattice. 

Evaluating weight multiplicities and Littlewood-Richardson coefficients is a 
difficult task. For At, computing Kostka numbers is immediate and Clebsch- 
Gordan's formula gives Littlewood-Richardson coefficients. For A2, one can still 
compute some small examples. But for general X r (r > 3) or for weights which 
components are big (say, with two digits), direct computation is usually intractable. 

There exist many formulae from representation theory for and Cyf^. The 
first one, valid in any complex semi-simple Lie algebra g, is Weyl's character 
formula 

ch(F(A)) = ^t£, where A^ = ^ (_i)=(-») e »0*) ) 
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where p is half the sum of positive roots for q. Littlewood-Richardson coefficients 
are obtained from this formula, since the character of V(X) ® V(p) is 



eh(V(\)®V((x)) = ch(V(A)) x ch(V(//)) = £ c x ^ch{V{v)). 

veL ; A+fi-i;eZ[A + ] 

But these two formulae do not lead to efficient computations when the rank of g or 
the size of coefficients of weights grow. Moreover, computing the whole character 
is untractable: for q — A 3 — s[4(C) and A = (2, 1,0, —3), the character ch(V(X)) 
has 9 monomials but the character ch(V(10A)) has 2903 monomials. 

Let us describe Kostant's and Steinberg's formulae in the case of any semi- 
simple Lie algebra q. Denote by k g (a) the number of ways one can write a vector 
a as a nonnegative linear combination of positive roots. Remark that k g (a) = 
unless a is in the root lattice Z[A + ]. This number satisfies the equation 



IL s a+(1 e «) a£z[A+] 

Let A and p be respectively a dominant weight and a weight such that A — p G 
Z[A + ]. A Weyl group element w <E W is valid for A and p if the root lattice 
element w(X + p) — (p + p) is in the cone C(A + ). The set of such u>'s is denoted by 
Val(A, p). Then Kostant's formula asserts that the weight multiplicity equals 



(2.1) c»= Yl (-lY {w) k s (w(\ + p)-(p + p)). 

tueVal(A,/j) 

Similarly let A, p, v, be three dominant weights such that A + p — v G Z[A + ]. 
The couple (w,w r ) £ W x W is valid for A, p, v, if the root lattice element 
w(X + p) + w'{p + p) — {y + 2p) is in C(A+). The set of such couples is denoted 
by Val(A,/i, v). Then Steinberg's formula asserts that the Littlewood-Richardson 
coefficient equals 



(2.2) c x \= (-iy {w)+£{w,) k B (w(\ + p)+w'(p + p)-(v + 2p)). 

( , w;,tu / )GVal(A,jU,y) 

Sets of valid Weyl group elements and valid couples of Weyl group elements turn 
out to be relatively small, when compared to W and W x W (which size is expo- 
nential in the rank). Remark that Kostant's (resp. Steinberg's) formula also work 
when A — p (resp. A + p — v) is not in the root lattice, since Kostant partition 
function vanishes on vectors that are not in the root lattice. 

From now on, let X r be a classical Lie algebra of rank r. Here X stands for 
A, B, C, D. Its positive roots system will be denoted by X+. 

Multiplicities and c Xjl behave nicely, in function of the parameters. More 
precisely, there exists a decomposition of the space t* © t* © t* in union of closed 
cones C, such that the restriction of c x v to each cone C is given by a quasi- 
polynomial function. This follows from theorems of Knutson-Tao KT99 (for 
A r ), Berenstein-Zelevinsky jBZOl) (for any semi-simple Lie algebra) giving 
as the number of points in a rational convex polytope. In the case of A r , the 
fact tha t c-y^ is given on each cone C by a polynomial function is proven in 
Rassart Rasfll], and the case of A 3 is treated as an illustration. The description 
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of the decomposition of t* © t* in cones C , where the function is polynomial 
for A r , was given for low ranks by Billey-Guillemin-Rassart BGR03 . See also 
Rassart's website [H] for wonderful slides. 

The common point to Kostant's and Steinberg's formulae is the function count- 
ing the number of decompositions of a root lattice element as a linear combination 
with nonnegative integral coefficients of positive roots of the Lie algebra. The next 
section deals with an efficient method to compute it. 

3. Counting integral points in rational convex polytopes 

3.1. Vector partition function. Let E ~ W and <E> be an integral matrix 
with set of columns A + = {(pi, . . . ,^>n} C E*. Choose a G Z' r . The rational 
convex polyhedron associated to $ and a is 

P($,a) = jx Y^Xifc = a, x t > oj . 

Remark 3.1. Every convex polyhedron can be realized under the form P($, a), 
that is as a set satisfying equality constraints on nonnegative variables. Indeed any 
inequality can be replaced by an equality by adding a new variable. For example 
polytopes {(x, y) G M 2 ; x > 0, y > 0, x + y < 1} and {(x, y, z) G R 3 ; x > 0, y > 
0, z > 0, x + y + z = 1} are isomorphic and have the same number of integral 
points. 

We assume that a is in the cone C($) spanned by nonnegative linear com- 
binations of the vectors 4>i, so that P(Q,a) in non-empty. We also assume that 
the kernel of $ intersects trivially with the positive orthant , so that the cone 
C($) is acute and P($,a) is a polytope (i.e. bounded). Finally, we assume that 
$ has rank r. The vector partition function is by definition 



k($,a) = |P($,a) nZ^| , 

that is the number of nonnegative integral solutions (xi, . . . , xjv) of the equation 
Eili = If $ = §>(X r ) is the matrix which columns are positive roots for a 
classical Lie algebra X r , then a i— > k($>(X r ),a) is the Kostant partition function. 
For example 



$(A 2 ) = -10 1 and $(S 2 ) = 

V -1 -1 J 

Note that the matrix &(A r ) has rank r (and not r + 1), since sums on lines are 
zero. 

A basic subset of A + is a basis o~ — {ai,...,a r } of E* constituted with 
elements of A + . Let B(A + ) be the collection of all basic subsets of A + . For such 
a a, let C(o~) be the cone of linear combinations with nonnegative coefficients of 
en's. Denote by Sing(A + ) the reunion of the facets of cones C(a), a G B(A + ); 
this is the set of singular vectors. Let C rog (A + ) := C(A + ) \Sing(A + ) be the set of 
regular vectors. A combinatorial chamber c is by definition a connected component 
of C reg (A + ). Combinatorial chambers are regions of quasi-polynomiality of the 
vector partition function a i— > fc($,a). Figure ^ represents cones C(A^) and 
C(B^), and their chamber decompositions. 



1 1 1 

V -i 1 o i ) 
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3.2. Brion-Szenes-Vergne formula for classical Lie algebras. Let us 

describe the formula, computing the number of integral points in rational convex 
polytopes P($>(X r ), a) associated to a classical algebra X r , that was implemented 
in our program. 

Let E = t and consider the set A + of positive roots for X r . Denote by A the 
set A + U (— A + ) of all roots. Let Ra be the vector space of fractions with poles 
on the hyperplanes defined as kernels of forms a £ A. Let 5a be the vector space 
generated by fractions f a := jt—^ — , a € B(A + ). Brion-Vergne BV97 proved 

that Ra decomposes as the direct sum Sa®9(Ra)- We define the Jeffrey- Kirwan 
residue of the chamber c as the linear form JK C on Sa ■ 



JK C (/ CT ) := 



voI(ct) -1 , ifccC(o-), 



0, if c n C(o-) = 0, 

where vol(cr) is the volume of the parallelopiped J^aeai®' ^\ a - We extend the JK 
residue to a linear form on Ra by setting it to on <9(-Ra), and to a linear form on 
the space of formal series Ra by setting it to on homogeneous elements of degree 
different from — r. For example, for the system A + = {ei, e<2, e\ +e2, &\ — e^\ Cl 2 
of positive roots for B 2 and the chamber c = l J+ e 1 © Z + (ei + e?) we have 

JK C ( ^) = JK C f ^) = JK C ( 1 - = -1, 
\xy z J \ xy z J \y z xy J 

since c C C({ei, e 2 }). 

Let T be the torus E/E%, where E% C E is the dual of the root lattice. Given a 

basic subset a, we define T(er) as the set of elements g & T such that e ( Q ^ 2l7rG '> = i 

for all «£<;; here G is a representative of g £ E/E%. Now let 

e (a.2inG+u) 

H9,a)(u) : =n QeA(1 _ e -( a , W) - 

Theorem 3.2 (Brion-Szenes-Vergne |BV99L ISV04j ). Let F C T be a finite 
set such that T{a) C F for all a € B(A + ). Fix a combinatorial chamber c. Then 
for all a £ Z[A + ] fie, we have: 

k(<l>,a) = Y,JK c (T(g,a)). 

gdF 

Now that we linked vector partition function and Jeffrey-Kirwan residue, we 
describe in Section 0] an efficient way to compute the latter. 
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4. DeConcini-Procesi 's maximal nested sets (MNS) DCP04 

We keep the same notations as in Section A subset S C A + is complete 
if S = (S) (~l A + . A complete subset is reducible if one can find a decomposition 
E = Ei © E 2 such that S = Si U S 2 with Si C E x and S 2 C E 2 ; else S is said 
irreducible. Let T be the collection of irreducible subsets. 

A collection M = {Ii,I 2 , ...,/„} of irreducible subsets Ij of A + is nested, if: 
for every subset {5*1, . . . , S m } of M such that there exist no i, j with Si C Sj, the 
union Si U S2 U • • • U S rn is complete and the Si's are its irreducible components. 
Note that a maximal nested set (MNS in short) has exactly r elements. 

Assume A + irreductible and fix a total order on it. For M = {Ii, . . . ,I S }, 
Ij G A + , take for every j the maximal element (3j G Ij . This defines an application 
4>{M) :— . . . , (3 S } C A + . A maximal nested set M is proper if <fi(M) is a basis 
of E*. Denote by V the collection of maximal proper nested sets (MPNS in short). 
We sort 4>(M) and get an ordered list 6(M) = [cti, . . . , a r ] . Thus 9 is an application 
from the collection of MPNSs to the collection of ordered basis of E* . For a given 
M, let then 



C(M) := C(ai,...,a r ), 
vol(M) := vol(e[ =1 [0,lK-), 
IResAf := Res Qr= o ■ • ■ Res Ql= o . 

Example 4.1. Let be the canonical basis of M. r , with dual basis e 1 (i = 
1, . . . , r), and define E as the subspace of vectors which sum of coordinates 
vanish. Consider the set A + = {e l — e 3 \ 1 < i < j < r} of positive roots for 
A r _i. Irreducible subsets of A + are indexed by subsets S of {1,2, ... ,r}, the 
corresponding irreducible subset being {e l — e J \i,j G S,i < j}. For instance 
S = {1, 2, 4} parametrizes the set of roots given by {e 1 — e 2 , e 2 — e 4 , e 1 — e 4 }. 

A nested set is represented by a collection M — {Si, S 2 , . . . ,Sk} of subsets of 
{1,2, ... ,r} such that if Si, Sj G M then either Si fl Sj is empty, or one of them 
is contained in another. 

For example one can easily compute that for the set of positive roots for A3 
(see Figure ^) there are only 7 MPNS, namely 



Mi = {[1, 2], [1, 2, 3], [1, 2, 3, 4]}, M 2 = {[2, 3], [1, 2, 3], [1, 2, 3, 4]}, 
M 3 = {[2, 3], [2, 3, 4], [1,2, 3, 4]}, M 4 = {[3,4], [2,3,4], [1,2,3,4]}, 
M 5 = {[1,3], [2, 4], [1, 2, 3, 4]}, M 6 = {[1, 2], [3, 4], [1, 2, 3, 4]}. 

Now we can quote the Theorem for the Jeffrey-Kirwan residue computation: 

Theorem 4.2 (DeConcini-Procesi). Let c be a combinatorial chamber and fix 
f G i?A- Take any regular vector u£t. Then: 

JK c (/)= £ ^) IRes M(/). 

Mev -. veC{M) y ' 

See [BBCV05] for a detailed description of how formulae from Theorems 13.21 
and 14. 21 were implemented. 

5. Our programs 

5.1. Description and implementation. Initial data for weight multiplic- 
ity and Littlewood-Richardson coefficients are only vectors (respectively two and 
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three). Our programs work with weights represented in the canonical basis of E* , 
and not in the fundamental weights basis for X r . Translation between these two 
bases is performed via straightforward procedures FromFundaToCanoX(r ,v' ) and 
FromCanoToFundaX(r , v) (where one replaces X by A, B, C, D, according to the 
algebra) . 

Computation of the weight multiplicity and of the Littlewood-Richardson 
coefficient c A t M is done by typing in 

Mult iplicityX( lambda, mu) ; 
TensorProductX( lambda, mu,mi) ; 

where A, fi, v are suitable weights. The syntax for computing quasipolynomials 
is slightly different. Assume that we want to evaluate (A',/i') i— > c y in a neigh- 
borhood of a couple (A,//), and (X',fjf, v') i— > c w v , in a neighborhood of a triple 
(A,/x,z/). Let Xp = [afi,...,a; r ], = [yi,...,y r ], v F = [z\,...,z r ], be three 
formal vectors where x^s, y^s and Zj's are variables. Then we use the command 
lines 

PolynomialMultiplicityX(lambda, lambdaF, mu, muF); 
PolynomialTensorProductX(lambda, lambdaF, mu, muF, nu, nuF); 

So for the polynomial (A', //) i-> c{, with A = (3, 2, 1, -6) and fi = (2, 2, -2, -2) 
for A3 we enter 

PolynomialMultiplicityA( 

[3, 2, 1, -6], [x[l], x[2], x[3], x[4]], [2, 2, -2, -2], [y[l], y[2], y[3], y[4]]); 

and get instantly 

i(3zi - 2yi + l)(3a:i - 2 yi + 2)(3s 1 + 6x 2 - 2y : - 6y 2 + 3). 


Remark that quasipolynomials c'^ and c t ^ v t ^ are obtained by setting X4 = t\i, 
yi = t/j,i, Z{ = tvi, so that 

PolynomialMultiplicityA( 

[3, 2, 1, -6], [3t, 2t, t, -6t], [2, 2, -2, -2], [2t, 2t, -2t, -2t]); 

returns (i + l)(t + 2)(t + 3)/6. 

Now some words about implementation. There are two main parts in our 
programs. The first one is the implementation of Theorems 13.21 and 14.21 it is 
described in [BBCV05]. The second one is the implementation of Kostant's l|2.1(l 
and Steinberg's (|2.1|l formulae using valid Weyl group elements and valid couples 
of Weyl group elements; it is a generalization for classical Lie algebras of what has 
been done for A r in |C03j . 

5.2. Comparative tests. Figure El describes efficiency area of the software 
LE and of our programs using MNS; any area located to the left of a colored 
line represents the range where a program can compute examples in a reasonable 
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time. Figures|3Hlpresent precise comparative tests of the software LIE, of DeLoera- 
McAllister's script |DM05j using LattE jLj and of our programs using MNS. 
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Figure 2. To the left, comparison for tensor product coefficients 
for A r : with LiE, with Sp(a) and with MNS. To the right, com- 
parison for weight multiplicity of a weight for B r : with LIE and 
with MNS. Similar Figures for C r and D r . 

All examples were runned on the same computer, a Pentium IV 1,13GHz with 
2Go of RAM memory. Remark that computation times for LattE and LIE are 
slower than those shown in DM05 , due to different computers. However, we 
performed exactly same examples for comparison purposes. 

As in DM05 , in Tables HUB] weights arc for jj[ r+1 (C) and not s[ r +i(C) (co- 
ordinates do not add to zero). However the sum of coordinates of A + fi — v 
vanish. 

Now some words about quasipolynomials computation. Let us examine the 
first example for B 3 in |DM05| . that is the evaluation of the quasipolynomial 
c t\ v tii f° r weights A = [0,15,5], [i = [12,15,3] and v — [6,15,6] expressed in 
the basis of fundamental weights. In canonical basis, these data become A = 
(35/2,35/2,5/2), fx = (57/2,33/2,3/2), v = (24,18,3). The program using the 
MNS algorithm returns the quasipolynomial 
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in 1099,4s. On the other hand, the computation of the full quasipolynomial c^" 
with formal vectors [xi, X2, Xa], [2/1,2/2,2/3], [^1,^2,^3] leads to a 87 pages result, 
obtained in only 1158,6s. With LattE, on our computer, one obtains the quasipoly- 
nomial c tA *" tM in only 825,8s. 

As announced in the introduction, our program is really efficient for weights 
with huge coefficients. Note that in the particular case of A r the MNS algorithm 
allows us to compute examples one rank further than the Sp(a) algorithm. 

The translation of the program using MNS in the language of the symbolic 
calculation software MuPAD is in progress. A version using distributed calculation 
on a grid of computers is in the air; it will considerably increase the speed of 
computations. 
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Figure 3. For A r , comparison of running times between the MNS algorithm, LattE and LiE 
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Figure 4. For A r , comparison of running times for large weights between the MNS algorithm and LattE 
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